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ABSTRACT 

I employ an ensemble of hydrodynamical simulations and the XSPEC MEKAL emission 
model to reproduce observable spectral and flux-weighted temperatures for 24 clusters. 
Each cluster is imaged at 16 points in its history, which allows the investigation of 
evolutionary effects on the mass-temperature relation. In the zero redshift scaling rela- 
tions, I find no evidence for a relationship between cluster temperature and formation 
epoch for those clusters which acquired 75% of their final mass since a redshift of 0.6. 
This result holds for both observable and intrinsic intracluster medium temperatures, 
and implies that halo formation epochs are not an important variable in analysis of 
observable cluster temperature functions. 
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1 INTRODUCTION 

Galaxy clusters are the youngest and largest organized struc- 
tures known to exist, and their ensemble properties can shed 
light on many cosmological problems. They arise from signif- 
icantly overdense regions on cosmological scales, which are 
exponentially rare events in a Gaussian initial perturbation 
spectrum. The shape and normalisation of the cluster mass 
function are therefore extremely sensitive to the statistical 
properties of the primordial density field. The evolution of 
cluster number densities is also tightly coupled to the cosmic 
expansion rate, and can precisely constrain fl m . Much effort 
has been expended towards a measurement of the cluster 
mass function and its evolution (see Henry 2000 for a sum- 
mary of recent work in this area), but much remains to be 
understood about the clusters themselves before we can be 
confident of such results. 

The properties of the local mass function can be con- 
strained by using either standard Press-Schecter (1974) (also 
Bond et al. 1991, Lacey & Cole 1993) theory, a more so- 
phisticated analytical model of the cosmic mass function 
(e.g. Sheth, Mo, & Tormen 1999), or fitting formula based 
on large-scale structure simulations (Jenkins et al. 2000). 
These methods predict the number density of dark matter 
haloes as a function of mass and redshift. Some relationship 
between model variables (i.e., the total mass within some 
density threshold) and a more easily observed quantity (e.g. 
X-ray temperature, X-ray luminosity, or weak lensing mass) 
must therefore be assumed to match cosmological predic- 



tions to the results of cluster surveys. Intracluster medium 
(ICM) X-ray temperatures show particular promise in this 
regard, since they demonstrate a tight (< 20% scatter) cor- 
relation with cluster mass components in both simulations 
(Mathiesen & Evrard 2001, total mass) and observations 
(Mohr, Mathiesen, & Evrard 1999, ICM mass), and the 
cluster X-ray temperature function (XTF) currently pro- 
vides the tightest constraints on the cluster mass function. 
It therefore seems likely that our best structure formation 
measurements of tt m will rely on medium-redshift measure- 
ments of the XTF, and it is essential to make sure that our 
interpretation of cluster temperature functions is correct. 

A modified form of the mass function allowing for dif- 
ferences between a cluster's formation redshift and observed 
redshift was first proposed by Kitayama & Suto (1996), and 
has begun to be commonly implemented in deriving con- 
straints on the power spectrum normalization erg at low red- 
shfits (Kitayama & Suto 1997, Kay & Bower 1999, Viana 
& Liddle 1999). This extension to the theory produces lit- 
tle change in the shape of the predicted mass function, but 
can have an appreciable effect on the predicted tempera- 
ture function if one assumes that clusters scatter around a 
mass-temperature scaling relation appropriate to their for- 
mation redshift. This implies that clusters have objectively 
identifiable formation events, during which the majority of 
a cluster's mass coalesces for the first time and the ICM 
virialises at a temperature appropriate to the halo mass and 
formation epoch. These assumptions are reasonable in the 
framework of standard Press-Schechter theory, but perhaps 
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pay too little heed to the more complicated, hydrodynamical 
evolution of the ICM. 

The goal of this paper is to investigate the evolution 
of intrinsic and observable cluster temperatures using an 
ensemble of 24 hydrodynamic cluster simulations, and test 
the appropriateness of this extension to the interpretation 
of cluster temperature functions. In the following section I 
describe the ensemble of simulated clusters and our model 
for the ICM X-ray emission. In section 3, I look for evi- 
dence of temperature evolution in our ensemble. Section 4 
summarizes the results of this letter. The Hubble constant 
is parameterized as Ho = lOO/i km s _1 Mpc -1 . 



2 SIMULATED AND OBSERVABLE 
TEMPERATURES 

We use an ensemble of 24 hydrodynamical cluster simula- 
tions, divided evenly between two reasonable cold dark mat- 
ter (CDM) cosmological models. These models are OCDM 
(fi = 0.3, cr 8 = 1.0, h = 0.8, T = 0.24) and ACDM 
(fio = 0.3, A = 0.7, a 8 = 1.0, h = 0.8, P = 0.24). Here cr 8 is 
the linearly evolved, present day power spectrum normaliza- 
tion on 8ft -1 Mpc scales. The initial conditions are Gaussian 
random fields consistent with a CDM transfer function with 
the specified V (e.g. Bond & Efstathiou 1984). The baryon 
density is set in each case to 20% of the total mass density 
(Q b = 0.2fi ). The simulation scheme is P3MSPH: the first 
stage is a P 3 M (dark matter only) simulation to find clus- 
ter formation sites in a large volume, and the second stage 
resimulates the formation of individual clusters with higher 
resolution. Smoothed particle hydrodynamics (SPH) is in- 
cluded in the individual cluster simulations to resolve the 
ICM structure in detail. The baryonic component is mod- 
eled with 32 3 particles, providing a typical mass resolution 
of 0.01% within the virial radius. The resulting sample cov- 
ers a little more than a decade in mass, ranging from about 
10 14 to 3 x 1O 15 M . These simulations were first presented a 
paper describing the cluster size-temperature relation (Mohr 
& Evrard 1997). 

The simulations model the dynamical and thermody- 
namical effects of gravitation, shock heating and adiabatic 
work on the ICM. Several potentially important pieces of 
physics are neglected. Radiative cooling is perhaps the most 
significant; our clusters cannot produce cooling flows in their 
cores. Cooling flows have the potential to greatly influence 
ICM luminosity and temperature measurements, but the en- 
ergy lost due to radiation is small compared to that released 
in the process of gravitational collapse. We therefore expect 
that the results of these simulations are comparable to obser- 
vational results which attempt to empirically account for the 
presence of cooling flows, either through excision or explicit 
modeling of excess core emission. Other neglected processes 
include galaxy feedback (Metzler & Evrard 1994) which can 
produce abundance gradients and shallower gas profiles; pre- 
heating of the ICM (Cavaliere, Menci, & Tozzi 1998; Lloyd- 
Davies, Ponman, & Cannon 2000), which can raise the ICM 
entropy and limit the density of baryonic cores; and elec- 
tron temperature lag, which slightly cools X-ray spectra in 



rich clusters (Chieze, Alimi & Teyssier 1998; Takizawa 1999) 
relative to the ion temperature. Further discussion of these 
issues in terms of their relevance to cluster simulations can 
be found in Mathiesen & Evrard (2001, hereafter ME01) and 
Bialek, Evrard, & Mohr (2001). 

In ME01 we introduced an ensemble of spectrally and 
spatially resolved X-ray surface brightness images derived 
from these simulations. We used the MEKAL (Mewe, Lemen, 
& van den Oord 1986) program from the XSPEC utility, 
since this is the emission model most commonly used to 
interpret observed ICM spectra. Each SPH particle was as- 
signed a 0.3 solar metallicity spectrum scaled to its density 
and thermodynamic temperature and binned in 50 eV inter- 
vals over the [0.1,20] keV bandpass. The clusters were then 
"observed" by collecting photons in a circular window cen- 
tered on the minimum potential of the cluster, producing 
a combined spectrum which incorporates emission from gas 
with a broad distribution of phases. We produced Chandra- 
like combined spectra and spectral images by adopting 150 
eV bins, a 0.5-9.5 keV bandpass, and finally convolving the 
photon distribution with Chandra's effective area function 
and a moderate (3.4 x 10 20 cm -2 ) absorbing column den- 
sity. The physical scale of the observation windows used in 
this paper varies from cluster to cluster, but usually corre- 
sponds to a fixed mean interior overdensity of 500 times the 
critical density appropriate to the redshift and cosmological 
model. This radius is labelled rsoo, and is a fixed fraction of 
the virial radius in any cosmology. The radius V2oo is also 
used in this work, but corresponds less closely to observable 
regions of the ICM. 

The resulting spectra are surprisingly similar in charac- 
ter to isothermal spectra with a temperature typically 10- 
20% lower than the mass-weighted mean thermodynamic 
ICM temperature. A semantic separation of the two con- 
cepts appears to be necessary, so hereafter we refer to the 
Chandra- like temperature just described as a spectral tem- 
perature, T s . The mass- weighted mean thermodynamic tem- 
perature will be referred to as the virial temperature T„, 
because this measure is found to precisely follow the virial 
relation Mtot oc T 1 ' 5 in these and other simulations. T v is 
identical to the mass-weighted temperature T m defined in 
ME01. When spectral fitting is limited to the 2.0-9.5 keV 
band (similar to most published temperature determina- 
tions) the deviation ST S between spectral and virial tem- 
peratures follows the relation 

ST S = ~ = (0.19±0.06) log 10 T s [keV]-(0.02±0.04).(l) 

This deviation steepens the observable mass-temperature re- 
lation to Aftot oc T 1 ' 62 , and implies that rich clusters are 
more massive than their spectral temperatures would lead 
us to believe. The difference between spectral and virial tem- 
peratures arises partly from the presence of recent minor 
mergers in the gas column, and partly from global ICM 
temperature profiles. Both effects contribute cooler gas to 
the observation window and an overabundance of soft, line- 
emission photons to the spectrum. I refer the reader to ME01 
for further details on this work. 

The temperature bias described in Equation hi has a 
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potentially important effect on any interpretation of the 
cluster temperature function and its evolution. As has al- 
ready been stated, these simulations do not include radiative 
cooling and are therefore directly comparable only to X-ray 
data which has accounted for the presence of cooling flows 
through excision or explicit modeling of the excess emission. 
Since it has not been possible until very recently to measure 
the spatial extent of cooling flows in high-redshift clusters, 
most studies of cluster evolution have evaluated the tem- 
perature and luminosity functions at low redshifts without 
attempting such corrections. In such works, the additional 
scatter introduced into cluster scaling relations by cooling 
flows has been accepted as a source of uncertainty in the 
cosmological constraints. 

With the advent of Chandra and XMM, however, we 
can do better. Markevitch (1998) has shown that by excis- 
ing the central 50/i _1 Mpc of nearby clusters and including 
a cooling flow component in the core spectrum, the scatter 
in the luminosity-temperature relation can be reduced by a 
factor of two. This allows for a more robust calculation of 
the maximum observable volume for each cluster as a func- 
tion of temperature, as well as providing a cleaner estimate 
of the ICM temperature in cooling flow clusters. The new 
telescopes should allow similar corrections to be made in 
high-redshift clusters. In order to achieve precise cosmologi- 
cal constraints from a measurement of evolution in the XTF, 
it is therefore desirable to model such cooling flow-corrected 
temperatures. 

Temperatures derived from high-quality ASCA data 
(Markevitch et al. 1998) are commonly used in constructing 
the local XTF or measuring the slope of cluster scaling rela- 
tions. These are flux-weighted mean spectral temperatures, 
which are biased towards regions containing dense and/or 
cool gas. This measure was found by Markevitch to produce 
a significant shift in cooling flow cluster temperatures rela- 
tive to an isothermal fit to the combined spectrum. Other 
cluster temperatures, on the other hand, were unchanged. 

Such temperatures are difficult to reproduce precisely 
in the simulations, mainly because the spatial extent and 
subdivision of ASCA images is different for each cluster in 
the sample (Markevitch et al. 1998) The large scale of ASCA 
spectral regions implies that the temperatures in each pixel 
are similar in character to the spectral temperatures de- 
scribed earlier, but with more weight given to the luminous 
core than the outer regions. I simulate these temperatures 
by dividing our observation windows into nine sectors with 
a morphology typical for ASCA clusters: a core region with 
radius rsoo /4; an inner annulus surrounding this region with 
outer radius 5rsoo/8, and an outer annulus surrounding this 
region with outer radius r 5 oo- The two annuli are each di- 
vided into 4 quadrants. The average spectral temperature T s 
over all nine sectors, weighted by total energy flux, will be 
referred to as a flux-weighted temperature Tf for the remain- 
der of this paper. This definition of Tf should not be con- 
fused with the emission-weighted temperature T e described 
in ME01, which is the density- weighted average thermody- 
namic temperature calculated over a spherical volume. Tf is, 
however, closely comparable to the emission-weighted tem- 
peratures described in the Markevitch et al. (1998) analysis. 



The deviation STf has a slope similar to that of ST S 
(Equation [lj), but a different normalisation: 

ST f = Tv ~ Tf = (0.22±0.05) log 10 2>[keV]-(0.11±0.03).(2) 

Tf is generally higher than T„ because of the extra weight 
given to the core, but the scale-dependence of its deviation is 
similar. The mean value of STf ranges from about -10% for 
poor clusters to +10% for rich clusters with virial tempera- 
tures on the order of 10 keV. Calculating the flux-weighted 
mean temperature for a coarse grid does not free tempera- 
ture measurements from a spectral bias due to the presence 
of multiple phases in a gas column, but it provides a more 
accurate estimate of the virial temperature for rich clusters. 

The deviations between virial temperatures and observ- 
able spectral temperatures described in Equations |l| and ^ 
do not vary with cosmological model or evolutionary epoch, 
nor do they vary significantly when the radius of the ob- 
served region is increased to enclose an mean overdensity of 
200 times the critical density. ST a and STf are plotted agains 
T 3 and Tf in Figure [jj These are robust measures of an ob- 
servational bias arising from realistic density and tempera- 
ture distributions in the ICM. SPH simulations have been 
shown to accurately reproduce the large-scale morphology 
of real clusters (Mohr et al. 1995), and a similar analysis of 
this particular ensemble displayed an even closer structural 
correspondence (Mohr & Evrard, private communication). 
N-body simulations of dark matter evolution have likewise 
been shown to produce merger histories which are in good 
agreement with Press-Schechter theory (Lacey & Cole 1994), 
so it is likely that these variations are similar in magnitude 
to those in real clusters. Analysis of Eulerian simulations 
reveals a similar level of clumping in the ICM (Bryan & 
Norman, private communication). 



3 MODELING THE MASS FUNCTION 

When converting a model mass function n(M, z) to an ob- 
servable XTF n(T, z), it is generally assumed that the clus- 
ters formed at their observed redshift. Kitayama & Suto 
(1996) were the first relax this assumption in modeling 
the XTF, applying the techniques of Lacey & Cole (1993) 
to calculate realistic distributions of formation epochs for 
the clusters at a given redshift. They reasoned that clus- 
ters which formed significantly earlier than their observation 
epoch would be hotter than clusters of a similar mass which 
collapsed more recently. Cosmological scaling of the back- 
ground mean density and temperature predicts a normal- 
isation evolution in the mass-temperature relation propor- 
tional to h(z)/ho = yj Sl m (l + z) 3 + for a flat universe 
(e.g. Bryan & Norman 1998), and the temperature correc- 
tions for such clusters can therefore be rather large. Other 
groups have recently begun to take up this standard in mod- 
eling the temperature function and its evolution (Viana & 
Liddle 1999, Kay & Bower 1999), but this extension to the 
model is not yet justified by observations. 

Making this correction implies an additional assump- 
tion, however: that the X-ray luminous regions of clusters 
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Figure 1. The upper panel displays the deviation 5T S between 
virial temperature and spectral temperature as a function of the 
spectral temperature. 384 points are displayed, corresponding to 
16 evolutionary epochs of 24 independent cluster simulations. The 
best-fitting correlation between 8T 3 and T s (Equation 1) is inde- 
pendent of redshift. The lower panel displays the corresponding 
relationship between 5Tj and Tt, as reported in Equation 2. The 
dashed lines in each panel display one-sigma variations on the 
best-fitting correlations. 

are approximately relaxed. Merger events have the potential 
to significantly alter a cluster's temperature when they oc- 
cur, and they need not be very large to do so (Cavaliere et 
al. 1999, ME01). For a cluster which formed at high redshift 
to maintain a temperature appropriate for that epoch, it 
should have already accumulated most of its observed mass. 
The rate of mergers observed in simulations makes this sce- 
nario seem unlikely, although Kitayama and Suto cite results 
from Eulerian simulations (Bryan & Norman 1998) indicat- 
ing that a cluster's temperature doesn't change much after 
its formation. We note that the temperature which they refer 
to is a simulation's luminosity-weighted temperature (sim- 
ilar to the T e used in ME01), which is very similar to the 
core temperature of the gas. Observational measures of the 
temperature such as T 3 and Tf are more heavily influenced 
by cool gas in the outer regions of the cluster, and more 
susceptible to minor merger events. 

These simulations can be used to test the sensitivity of 
observable and virial temperatures to ongoing minor merg- 
ers. The mass of our clusters as a function of flux-weighted 
temperature is plotted in Figure [| and presents a tight cor- 
relation with only 18% scatter: 

log 10 (Mtot/i(») = (1.66 ±0.04) log 1Q T> + (13.59 ±0.02).(3) 

This plot combines cluster outputs at redshifts z = TO, 0.5, 
and 0. The small degree of scatter in this plot strongly im- 
plies that there is no significant contamination of the tem- 
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Figure 2. The relationship between flux-weighted temperature 
Tf and total mass measured within rsoo- Clusters at redshifts 
0.0, 0.5 and 1.0 are included in the plot, with masses scaled by 
the cosmological evolution factor h(z). The best- fitting relation is 
drawn, and its parameters are reported in Equation H. 



perature ensemble by clusters which virialized early; when 
scaled to a similar epoch the three relations M(Tf,z) are 
identical. We can, however, probe this issue more deeply. 

We define a cluster's formation epoch Zf as the red- 
shift at which it has acquired 75% ± 7.5% of its final mass, 
following the convention used by Viana & Liddle (1999) in 
their recent paper constraining cosmological parameters. If 
the extension of Kitayama & Suto is relevant to our inter- 
pretation of the local XTF, then we should see a correlation 
between cluster temperature and Zf. clusters with high for- 
mation redshifts should, on the average, have temperatures 
higher than the mean mass-temperature relation. Figure ^ 
plots this difference for virial, spectral, and flux-weighted 
temperatures within observation windows of radius rsoo and 
r2oo- Error bars along the Zf axis are given to objects which 
passed through the mass threshold with a significant compo- 
nent of continuous accretion, so that the cluster had between 
67.5% and 82.5% of its final mass in more than one output 
frame. The uncertainties implied by these error bars are not 
used in the statistical analysis. All the best-fitting lines for 
the six data sets are consistent with zero slope, although it is 
fair to say that there is evidence for a slight positive correla- 
tion. On the other hand, this correlation is largely driven by 
the rare clusters which formed at very high redshifts; most of 
these objects formed at Zf < 0.6 and are evenly distributed 
about the mean mass-temperature relationship. 

Correlation coefficients and best-fit line parameters for 
these data are summarized in Table |l|. The correlations 
within r2oo are significantly influenced by a cluster with a 
formation redshift of 0.8 and very large error bars on that 
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Figure 3. The best-fitting mass-temperature relationship T(M) 
is calculated for the 24 clusters in this ensemble at a redshift 
of zero, and the deviation of individual clusters from this rela- 
tion AlogT = log(T/T) is plotted against their formation red- 
shift Zt . Three definitions of the temperature (virial, spectral, and 
flux-weighted) and two observation windows are examined. If the 
temperature of a cluster were strongly influenced by its forma- 
tion epoch, we would expect to see a correlation between A log T 
and Zf. The evidence for correlation in this data set is marginal 
at best, and negligible for the flux- weighted temperatures which 
are most similar to ASCA measurements. Although the simulated 
clusters display the expected cosmological evolution in the nor- 
malization of their mass-temperature relation, the scaling relation 
at zero redshift does not appear to be contaminated by clusters 
which virialized at earlier epochs. 



value. This cluster acquired most of its mass very early in the 
simulation and grew through gradual accretion thereafter, so 
the reshift range during which it had a mass of 75%±7.5% its 
final mass is very long. It's temperature, while higher than 
the ensemble average, is well within the variations seen for 
more recently formed clusters. If this point is left out the 
analysis, then all the slopes become consistent with zero in 
their one-sigma uncertainty range, and the correlation coef- 
ficients drop to 0.24 (T„), 0.19 (T s ), and 0.005 (T». These 
coefficients correspond respectively to 32%, 39%, and 100% 
probabilities of uncorrelated data. 

The discrepancy between our intuition (clusters which 
first virialized at an early epoch should be hotter) and these 
simulations can be resolved by acknowledging the essentially 
dynamic nature of clusters in a low-density universe. Multi- 
ple lines of observational evidence point to an fl m ~ 0.3 cos- 
mology in which clusters are still forming at the present day, 
and the theoretical construct of a relaxed, virialized cluster 
seems to have few counterparts in the observable popula- 
tion. Rather than treating clusters as static fossils of the 
primordial density field, we should attempt to model them 



Temperature 


R 


slope 


intercept 


^200 


T 


0.31(15%) 


0.060 ± 0.044 


-0.018 ±0.017 




T s 


0.28(20%) 


0.080 ± 0.066 


-0.024 ± 0.025 




T f 


0.12(58%) 


0.030 ± 0.060 


-0.009 ± 0.022 




r 


0.21(33%) 


0.041 ±0.051 


-0.014 ±0.020 




T a 


0.28(20%) 


0.078 ± 0.086 


-0.026 ± 0.034 






0.11(68%) 


0.024 ±0.057 


-0.008 ± 0.022 



Table 1. A correlation analysis of the data presented in Figure Q 
R is the correlation coefficient of the data, and is also translated 
into the probability that an uncorrelated set of 24 random points 
would produce a correlation coefficient at least that large. The 
slope and intercept of the best-fitting lines plotting in Figure ^ 
are also given. 



explicitly as evolving entities. One example of such a model 
has been presented by Cavaliere et al. (1999), who analyze 
ICM structure in terms of "punctuated equilibrium", a se- 
quence of merger shocks followed by partial relaxation of the 
ICM to the shock boundary conditions. Their model agrees 
with these simulations in showing that minor merger events 
should have a more important influence on the evolution 
of ICM temperatures than major mergers, and makes some 
important predictions about the behavior of ICM scaling 
relations for rich groups and poor clusters. 



4 CONCLUSIONS 

An analysis of the evolutionary history of simulated clusters 
shows only a slight dependence of cluster temperature on 
formation epoch out to z{ ~ 0.6, although the small num- 
ber of clusters in this ensemble leaves open the possibility 
of evolution in the temperatures of clusters which formed at 
high redshifts. It is also worth reiterating that when using a 
flux- weighted temperature Tf, which is the definition most 
appropriate to analysis of the cluster temperature function, 
this correlation essentially vanishes. This lack of dependence 
on formation epoch can be traced to a high frequency of 
smaller (mass ratio 25%) merger events, which introduce 
cool gas into the ICM and allow it to approach an equilib- 
rium appropriate to the merger epoch. This work implies 
that extending analysis of the temperature function to ac- 
count for the difference between a cluster's formation and 
observation redshifts is not necessary in a dynamically young 
halo population. 
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